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Abstract 



The R- and K-matrix parametrizations are analyzed and compared for the 
elastic a-a scattering at center-of-mass energies below 40 MeV. The two 
parametrizations differ in their definitions of the resonance energy which can 
lead to quite different results. The physical values of the best-fit parameters 
are compared with those computed for a potential model. The existence of a 
broad resonance near 9 MeV is not supported by the data or by the potential 
model. We discuss the positive and negative aspects for both parametriza- 
tions. 
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I. INTRODUCTION 



The astrophysical S-factor for the reaction ^^C(a,7)^^0 at the astrophysically most ef- 
fective energy of E=0.3 MeV has been obtained from the extrapolation of a parametrized 
cross section by fitting data in the center of mass energy range between 1 and 3 MeV. Al- 
though Azuma et al. have recently found R-matrix and K-matrix parametrizations which 
give nearly the same results this agreement has not always been observed and fits with 
quite distinct differences between K- and R-matrix have been found 0-0- This situation 
has motivated us to make detailed tests of the two parametrizations and to compare the 
physical values from the two parametrizations. To this end, one must turn to a much sim- 
pler problem than the simultaneous parametrization of the three data sets for ^^C(a;, 7)^^0, 
^^C(a,a)^^C and ^^N /3-decay. Following Barker [TU| we consider the s-wave a + a elastic 
scattering. This problem also has the advantage that, besides accurate data, an excellent 



potential model description [rT| is available which we can use as a benchmark. Below the 
energy of the first reaction threshold, ''Li+p at 17.3 MeV, the 43 data we use Qare the same 
as in reference fl^. But, when it is useful to consider higher energies, we supplement them 
by the real parts of 5 complex phase shifts obtained by Darriulat et al. |T2|, as has been 
done in in the energy range 26-40 MeV. In the present paper, all energies refer to the 



center-of-mass system. 

The phase shifts corresponding to the potential model are obtained from the radial wave 
function M(r, E) which solves the radial Schrodinger equation 



^ {Vn + Vc-E) 



u{r,E) = 



subject to the boundary condition 



u{r = 0,E) =0 . 
Here the nuclear and Coulomb potentials have the form [jll 

VAr(r) = Voexp (-6r^) ; 



(1) 



(2) 



(3) 



Voir) = 4eVf (/3r)/r , (4) 

respectively. The best fit to the data is obtained with b = 0.212 fm~^ and (3 = 0.75 fm~^. 
The depth Vq will be chosen in such a way that the R-matrix and K-matrix parametrizations 
have a pole at the energy of the ^Be ground state, i.e. at the energy Eg = 92.08 keV [|T^]. 



Kindly communicated by F.C. Barker. 



2 



II. R-MATRIX FITS 



In terms of the R- matrix |]TE[, the s-wave phase shifts are 

where a is the channel radius, B the real boundary condition constant, and 

0(a, i?) = arctan(— ) , (6) 

S(a,E) = P{a,E)(F„Fi + G„Gi) . (8) 

— 0, P, and S" are the hardsphere phase shift, the penetration factor and the shift function, 
respectively. The radial Coulomb wave functions Fq, Gq depend on 



P 



kr = 0.309428rFi/2; r] = 0.891132^-^/2 ^ (9) 



where r is in fm and E is in MeV, while in Eq. (8) the primes stand for d/dp. For the sake 
of comparison we adopt the radius parameter a = 6 fm from Ref. [|10|. (We note that this 



value appears to be somewhat large. For example. Fowler et al. used a = 4 fm and the 

1 /3 1 /3 

conventional value is a = 1.44(Ai + A2 ) = 4.57 fm. On the other hand, with a = 6 fm, 
the phase shifts have not quite reached their asymptotic values, as we will see below.) 
In terms of the radial factor u{r, E), the R-matrix is 

_ /^O) 

[a{du/ dr)r=a — Bu{a, E)] 

Let us call Ei,E2, ... the eigenenergies satisfying the boundary condition 

a{du/dr)r=a — Bu{a, E) = . (11) 

Together with (1) and (2) this defines a classical Sturm-Liouville problem |14| whose eigen- 
values are all real. The R-matrix has been defined in such a way that it has poles at the 
-B-dependent eigenenergies Ei. Its expansion then reads 

00 2 



where the 7? are the formal []T3[ reduced widths, in terms of which the formal widths are 

T, = 2^^P{a,E) . (13) 
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Note that the Ei are not resonance energies. However, in R-matrix theory a resonance is 
associated with each eigenenergy. Let us first consider the case of a one-pole approximation, 
R = -ff/{Ei -E). We then have 



arctan 



r./2 



E, + ^,{a,E)-E^ 



where the energy shift is 



Ai{a,E)=^l[B~S{a,E)] . 



The resonance energy Ei^j. is defined as the shifted Ei satisfying the equation 113 



Ei + AAa.E,,r)-E, 



■'i,r 







(14) 



(15) 



(16) 



Assuming that in the neighborhood of Ei^^ a hnear approximation of Aj(a, E) is satisfactory, 

(17) 



the so-called observed |T3[ reduced width and observed width are 



^^^^ [l+^KdS/dE)E=E., 



rO = 2(7,Wa,i?) 



and we have 



(5 + = arctan 



r°/2 

Ei^r — E 



For the general theory following from Eqs. (5) and (10) we have 
S{E) + 0(a, E) = arctan [P(a, E)Q{a, E)] 

with 

1 u{a,E) 



Q{a,E) 



(i?-i + B- S) a{du/dr)r=a - S{a, E)u{a, E) 



(18) 



(19) 



(20) 



(21) 



The function Q{a, E) is neither an R-function in the sense of Wigner| nor a meromorphic 
function of E. As a generalization of Eq. (16), the resonances are defined as the real energies 
at which 



R-^ + B-S = 



(22) 



and we will call them Ej''\j = 1, 2, .... They are pole energies of Q{a, E) and 



2 See Ref. |13] p. 277 and the references there to E.P. Wigner's original papers on the mathematical 
R-functions. 
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5{Ef) + (j){a,Ef) =90° {mod 180°) . (23) 
The pole term of Q corresponding to ■* is 



with 



{Ef - E) 



(24) 



(25) 



^E^ ^i?=s: 

The Ei^r and (7^)^ defined by Eqs. (16) and (17) are one-level approximations of the gen- 
erahzed quantities and (7°)^ defined by Eqs. (22) and (25). When a good potential 
model is available, Q is given by Eq. (21). However, generally this is not the case. Provided 
an R-matrix fit has been obtained to the data, one can obtain the best values for E^''^ and 
{liY by substituting the parametrized R-function into Eqs. (20) to (25). 

When fitting data it is often convenient to choose the constant B so that the quantity 
B — S{a, E) vanishes in Eq. (5) at one of the energies e\^\ say at E^\ When B = S{a, E^^)^ 
we will call the eigenenergies which satisfy Eq. (11) by e'i\ E2'\ .... 

Turning to our example case, elastic a + a scattering, we choose e[^^ = Eg and B = 
S{6 fm,Eg). Then the boundary condition (11) is satisfied at E = Ei'\ if the potential 
depth is Vo = -119.216148 MeV. Since 0(a = 6 im, Eg) is very small (equal to 0.00029°), 
we have to a very good approximation 

6iEg) = 90° ■,{d6/dE)E=E,>0 (26) 

and Eg practically coincides with a resonance energy as defined by the conventional defini- 
tion. 

In Fig. 1 the phase shifts obtained with the potential model are compared with the 
data. How do the poles in the best-fit R-matrix compare with those of the potential model? 
To answer this question we have computed the energies the reduced and the formal 
widths for z = l — 3,j = l — 3 and the observed widths for i = 1 — 3. The results are 
in Table 1, while we have summarized the best-fit R-matrix parametrizations to the data 
below 18 MeV in Table 2, assuming B = S{a = 6 fm, -E'j"'^) with j = 1 — 3. No widths can 

be attributed to E^ and E^ , which are below threshold. But it is easily verified that, 
using the one-pole approximation with Eq. (16), El = —0.395 MeV in Table 2 is shifted 

to Ei,r = 0.101 MeV, close to the resonance energy, while using Eq. (22), Ei is shifted 

(3) 

exactly to the resonance energy Eg. Similar results hold for ^ = -0.379 MeV. In table 
2, the three ^^e different. This might be surprising at first glance, as a transformation 
of the R-matrix parameters from one boundary condition to another should not change the 
quality of the fits. However, this is only true if all R-matrix parameters are allowed to vary 
p!7| , which is not the case here. 

We also used the results in Table 1 to calculate the R-matrix phase shifts with a 3-pole 
approximation for the case B = S{a = 6 fm,E[^^). In Fig. 2, these calculated phase shifts 
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are compared with the exact potential model phase shifts. The agreement is not good. 
Comparing Tables 1 and 2, one observes that the reduced widths 7I of the third poles in 
Table 2 (2.112, 1.892, and 1.821 MeV) are much larger than those in Table 1 (0.734, 0.769, 
and 0.770 MeV). Clearly, in the best fits (in Table 2), the large reduced widths at the third 
poles compensate the contributions from poles at much higher energies. The convergence 
of the sum of poles calculated for the potential model is very slow. Even adding the fourth 
and fifth poles (at ^ = 56.0 MeV and E^^ = 89.2 MeV with the reduced widths 0.688 
and 0.660 MeV, respectively) is not sufficient for a good fit, as seen in Fig. 2. 

Like Q{a,E), the observed resonance energies E^''^ and the corresponding observed re- 
duced widths (7°)^ are -B-independent, but they also depend strongly on the radius a. This 
is illustrated in Table 3 and Fig. 3 for the potential model. In Fig. 3, with different channel 
radii, we plotted the phase 6 + (f) corresponding to the data and to the potential model. The 
best fit of the potential model to the 48 data is obtained with a = 5.5 fm and = 84.1. 
For each radius the potential depth Vq has been modified to satisfy exactly e[^^ = Eg. In 
the range a = 7 — 8 fm, the asymptotic phase has been reached and the fit remains good 
(x^ < 100). The density of resonances increases strongly with the channel radius without 
greatly changing x^- 

The uncertainty regarding the "proper" choice of the channel radius is such that, in 
practice, it is often chosen to give the best fit to the data. 



III. K-MATRIX FITS 



For the potential model the conventional K-matrix reads 

^ uF^ - u'Fo ^ 



and the phase shifts are 



'''' 



5(E) = arctanK . (28) 



With the channel radius a = 6 fm and the potential depth Vq = —119.217576 MeV, the 
K-matrix has a pole at ei = Eg = 92.08 keV, while other real poles are at 62 = 2.808641 MeV 
and at 63 = 31.881031 MeV. But K cannot be parametrized as a sum of pole terms, because 
of the essential singularities of the radial Coulomb wave functions Fo,Gq. The K-matrix has 
an infinite number of complex poles converging to E = 0. 

To eliminate these singularities, a modified /C-matrix has been defined in Ref. It 



has no other singularities than isolated poles. It is obtained by separating the threshold 
factors of Fq, Gq which depend only on 77, and substituting into Gq a polynomial in 77"^ (i.e. 
in E) to the singular function 



^ The potential depths Vq are not the same in the R- and K-matrix parametrizations. However, 
because the energy Eg is very small we have |Fo| << |Go| and S{a,Eg) ^ (pGq/Gq), so that the 
two depths are nearly the same. 
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h{r]) = ^^{l+tr]) + ^^{l-t7])-\nr] , (29) 

where is the digamma function. The function /i(?7) has an essential singularity at rj = oo, 
i.e. at E' = 0. At real energies, over a finite range (say for E between ri and r2), the 
polynomial is chosen to fit the function h{r]) for E > and the function 

h~^{r]) = '4'{ir]) + — ln(z?7) (30) 

for < 0. This is easily achieved, at any desired approximation, using Chebyshev polyno- 
mials up to the appropriate degree in rf'^ . Let n be the maximum degree chosen for these 
polynomials. For a given n, one can then choose r\ and r2 so that the polynomial in r/"^ is 
exactly equal to h{r]) at E = ei and 62- With n = 95 and 

ri = -4.874644MeV ; ra = 44.591714MeV (31) 

the real poles ei, 62 are the same for K and /C, and the 63 pole is only very slightly shifted 
to 63 = 31.881027 MeV. 

The energy ei of the first pole satisfies the usual resonance conditions 

S{e,) = 90°(mod 180°) , (^)i?=e. > (32) 

and the reduced width gf is positive (see Table 4). For the two other poles, the energies 62, 
and 63 satisfy only the first condition. Their reduced widths are negative. They are echo 
poles and contribute only to the background of /C. 

The /C-matrix corresponding to the potential model is the meromorphic function 

( uFj - u'F„ \ 



where Fo,Go are the modified Coulomb wave functions [T^, while for the phase shifts, we 
have 

S{E) = arctan(p2/C) (34) 

with 

^ = (exp{27rr;} - 1) ^^^^ 

The very way we have defined the /C matrix has several consequences. Contrary to K, /C 
can be expanded in a series of pole terms 



^ In order to have dimensionless /C and reduced widths with an energy dimension, we substitute 
^-1-1/2 £qj, /j'+i/2 following the definition of p by Eq. (4.7) in Ref. [16|. 
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oo 2 

i=i (ej - tj) 

At positive real energies, the phase shifts defined by Eqs. (28) and (34) are practically equal, 
while, at complex energies, K, is different from K/p^. Like the polynomial in ?7~^, K, also 
depends on the choice of n,ri, and r2- 

With the potential model, the three real poles of K, below 40 MeV and their reduced 
widths are given in Table 4. With the corresponding 3-pole approximation of /C, one does 
not obtain a good fit of the phase shifts. A better approximation requires more real and/or 
complex poles. Like the expansion of the R-matrix, the expansion of /C converges only very 
slowly, and the existence of complex poles introduces further complications. Moreover, only 
the real poles are uniquely defined, in the sense that the energy of the complex poles and 
their residues depend on the choice of n, ri and r2- In particular, the density of complex 
poles near the real axis increases when a larger n is chosen. 

In the domain -4 MeV < ReE < 44 MeV, - 5 MeV < \m.E < 5 MeV with n = 45, 70, 95 
the number of pairs of complex conjugate poles is 7, 22, and 37, respectively. Many of these 
poles have very small complex residues and their contributions to the /C-matrix are negligible. 

In a /C-matrix fit to the data, when an additional background term is needed, it cannot 
be a real pole in the energy range of the data, since 6 = 90° (mod 180°) only at ei, 62, 63. But 
a real background pole at negative energy cannot be excluded a priori. In fact, the energy 
dependence of the /C-matrix background plotted in Fig. 2 of [l^ shows a decreasing positive 



background with a concave curvature. Since the energy dependence of an echo pole below 
threshold has also a concave curvature at positive energies, introducing such a pole into the 
parametrization seems to be the simplest way to get a good fit below 18 MeV. With the 
Ci and 62 poles and a 3-parameter background composed of a constant and an echo pole at 
—0.327 MeV, Barker [|I^ easily obtained a very good fit to the data. We checked that the /C 
matrix for the potential model has no such pole at negative energies and we must conclude 
that the three parameters of the background are only ad hoc parameters without physical 
meaning. 

To further confirm this conjecture, we made another best fit to the data up to 40 MeV 
using the three real poles and a pair of complex conjugate poles. With 64 = e/j + ze/ and 
S'l = 7r + ^7/) the background term has the form 

+ conj. = — t^^o , 2 ■ (37) 



2{e,-E) (en-Ey + ej 

A very good fit is obtained to the 48 data points up to 40 MeV (x^ = 27.0). The results are 
in Table 4 and Fig. 4. They confirm that in a /C matrix fit to data with a minimal number 
of parameters, the background parameters are unlikely to correspond to the complex poles 
computed for a phase-equivalent potential model. 



IV. CONCLUSIONS 

The main advantage of an R-matrix parametrization is that R can be expanded in terms 
of real pole terms because R is an i?-function in the sense of Wigner. The parameters 
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introduced in an R-matrix fit are the eigenenergies Ei and the formal reduced widths 7^ . 
However, this is a parametrization of S + (f) rather than of 6. One sees in Fig. 3 that at 
E >2 MeV (5 + is steadily increasing with E and one can obtain the energies, at which 

5(E) + 0(a,E) = 90°(mod 180°) , (38) 

defining the so-called observed R-matrix resonance energies. The corresponding observed 
reduced widths (7°)^ and the observed widths r° are obtained from Eqs. (16) and (17), or, 
if a 1-pole approximation does not apply, from Eqs. (22) and (25). 

Do all these so-called resonance energies E^^^ really correspond to physical resonances? 
How are they compared to the resonances defined by the usual conditions 

6{E) = 90°(mod 180°), 37; > ? (39) 

dE 

Let us first evaluate the time delay in the scattering process, which is 2v~^dS/dk [|r^, where 
V = hk/M (M being the reduced mass) and 

dS d d 

— = — [arctan(P(a, E)Q{a, E))] - — 0(a, E) . (40) 

With the potential model, aX E = e'i \ the two terms on the right hand side of Eq. (40) 
are 61 ■ 10^ and 97 ■ 10"^ fm, respectively, and the resonance conditions (39) are satisfied. 
ki E = Ef\ the same terms are 3.7 and 6.6 fm, respectively. Thus, d6/dk is negative and 
this corresponds to a time advance, excluding a physical resonance. A similar argument 
holds also at E = Ef^ when i > 2, but in Ref. 0, the third pole at 32.9 MeV was already 

considered background term. 

(2) 

Let us now see how the E{ R-matrix pole appears as a complex pole of the iS-matrix. 
Designating that pole by p{E), we obtained 

, , (-0.0535 + 20.0339) , , 

= — (e:^) — 

where E^ = 3.34 - ilQ.59 MeV is far off the real axis. With 5(3.34) = -0.921 + zO.389 and 
p(3.34) = —0.00204 — iO. 00322, one finds that p{E) does not add a resonant contribution to 
SiE). 

Under such conditions, we must conclude that, in the best fit, E2'^ and E^^^ are both 
background poles, despite the fact that 7I and 7I (like 7^) are positive. This has been 
confusing [|l^ and, like the strong dependence of the on the channel radius, it is a weak 
point of the R-matrix parametrization. This does not concern the good quality of fits one 
can obtain with this method, but the physical interpretation of the parameters of the pole 
terms. 

In a /C matrix parametrization, the approximate energies at which there are real poles of 
K. are directly suggested by the data, since they are the energies at which 6 = 90° (mod 180°). 
This is the main advantage of a /C-matrix parametrization. These energies are either res- 
onances or echo pole energies according to whether the phase shift data are increasing or 
decreasing at the energies concerned. In the present case, the echo poles at 2.8 MeV and 
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31.8 MeV arc only a part of the background terms of /C. At least one more background 
term is needed to obtain a good fit to the data. It can be a constant, a real echo pole below 
threshold or at higher energies than the data, or a pair of conjugate complex energy poles. 
The parameters of these poles, energies and residues, are no more than ad hoc parameters. 
They are not expected to be poles of the modified /C-matrix corresponding to a potential 
model. This is the weak point of a /C-matrix parametrization. 

We can summarize as follows the lessons we have learned from the present analysis of 
the a + a scattering using the R-matrix and /C-matrix methods for pratical fit procedures. 
In order to obtain a good fit to data, one should use as fiexible parametrizations as possible, 
without worrying about the relation of certain parameters (e.g. echo poles below threshold 
in the /C-matrix or R-matrix eigenenergies and residues) to physical quantities (e.g. reso- 
nance parameters). In certain cases such a connection cannot be made and some of the fit 
parameters are purely ad hoc. Nevertheless, one of the main purposes of R-matrix and /C- 
matrix fits is to give a reasonable and physically motivated basis to extrapolate data. If both 
methods are used carefully, they can be expected to give similarly good parametrizations, 
as in the case of ^^C(q;, 7)^^0. 
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FIGURES 



FIG. 1. Comparison of the experimental phase shifts (points) with the potential model (solid 
line). The dashed line corresponds to the best R-matrix fits in Table 2. Note that the 3 best fits 
are indistinguishable within the thickness of the line. 

FIG. 2. Comparison of the potential model phase shifts (solid line) with the computed 3-pole 
(long-dashed) and 5-pole (short-dashed) R-matrix approximations. Compared with the data, the 
values are 93.9 (potential model), 9202 (3-pole) and 3897 (5-pole), respectively. 

FIG. 3. Values ior 5 + cf) for the data and the potential model to illustrate the dependence of 
the R-matrix ' resonances on the channel radius. 

FIG. 4. /C-matrix fits to the data. The solid curve (fitted to the data up to 40 MeV) has been 
obtained with a pair of conjugate complex poles as background, and the dashed line (fitted to the 
data up to 20 MeV) with a background consisting of a constant and an echo pole below threshold. 
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TABLES 



TABLE L With the potential-model, three different computed parametrizations of the 
R-matrix are obtained with the boundary-condition constants B = S{Q,eI^^^) {i = 1,2,3), re- 
spectively. All the parameters are in MeV. 



i 


1 


2 


3 


4 


5 




0.09208 


9.764 


29.274 


56.006 


89.217 


ll 


0.228 


0.835 


0.734 


0.688 


0.660 




8.15E-6 


9.207 


14.500 


18.947 


23.022 




-0.343 


8.419 


28.173 








0.386 


1.004 


0.769 










10.195 


14.903 








-0.352 


8.397 


28.156 








0.389 


1.006 


0.770 






Li 




10.203 


14.904 








0.174 


1.000 


0.769 






n 


6.20E-6 


10.156 


14.901 







TABLE n. Best R-matrix fits below 18 MeV. The parameters in parentheses denote fixed input 
values. The fixed parameters for the second and third best fits are obtained using Eq. (21) with 
the parametrized R-matrix from the first best fit. All the parameters are in MeV. The three 



are not identical because not all parameters have been varied |17|. 





1 


2 


3 






(0.09208) 


9.787 


32.954 






(0.199) 


0.848 


2.112 






(7.099E-6) 


9.366 


44.359 


18.05 




-0.395 


(8.424) 


30.463 






0.343 


(1.012) 


1.892 








(10.281) 


38.162 


15.86 


eP 


-0.379 


8.406 


(30.032) 




ll 


0.336 


1.010 


(1.821) 








10.243 


(36.459) 


15.90 


{iff 


(0.156) 


(1.008) 


(1.820) 






(5.57E-6) 


(10.241) 


(36.441) 
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TABLE III. Nuclear potential depths, and second and third poles of Q{a,E), for different 
values of a (in fm). The depths Vq (in MeV) are chosen so that = Eg for each radius a. Vq 
and the E^^^ are in MeV. The x^-values have been calculated for 48 data points. 



_a Vo El El 

4.0 -122.8034 28.755 78.050 1464 

5.0 -119.5106 14.801 44.871 143 

5.5 -119.2750 11.018 35.186 84.1 

6.0 -119.2161 8.418 28.155 93.9 

7.0 -119.2003 5.279 18.980 100 

8.0 -119.1999 3.595 13.516 99.7 



TABLE IV. /C-matrix parameters computed from the potential-model and two best fits up to 
18 and 40 MeV, respectively. Except for the dimensionless constant in the first best fit, all the 
parameters are in MeV. 





Potential i 


model 


Best fit 


Best fit 


Range 


0.01-40 MeV 


0.01-18 MeV 


0.01- 


-40 MeV 


j 
1 


0.09208 


46.138 


(0.09208) (45.726) 


(0.09208) 


(45.726) 


2 


2.809 


-10.801 


2.819 -10.278 


2.809 


-10.924 


3 


31.881 


-6.364 


-0.324 -39.504 


31.865 


-5.236 


4 






const. =-0.656 


-1.026±i0.199 


-31.082±i307.53 


Ti 


5.62E- 


-6 


(5.57E-6) 


(5.57E-6) 


x' 


92.5 




18.2 




27.0 



14 




2 4 6 8 10 12 14 16 18 20 



E (MeV) 




2 4 6 8 10 12 14 16 18 20 

E (MeV) 
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E (MeV) 
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E (MeV) 



